Turbulence induced additional deceleration in relativistic 
shock wave propagation: implications for gamma-ray 
burst 
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Abstract The late afterglow of gamma-ray burst is be- 
lieved to be due to progressive deceleration of the for- 
ward shock wave driven by the gamma-ray burst ejecta 
propagating in the interstellar medium. We study the 
dynamic effect of interstellar turbulence on shock wave 
propagation. It is shown that the shock wave deceler- 
ates more quickly than previously assumed without the 
turbulence. As an observational consequence, an earlier 
jet break will appear in the light curve of the forward 
shock wave. The scatter of the jet-corrected energy re- 
lease for gamma-ray burst, inferred from the jet-break, 
may be partly due to the physical uncertainties in the 
turbulence/shock wave interaction. This uncertainties 
also exist in two shell collisions in the well-known inter- 
nal shock model proposed for gamma-ray burst prompt 
emission. The large scatters of known luminosity re- 
lations of gamma-ray burst may be intrinsic and thus 
gamma-ray burst is not a good standard candle. We 
also discuss the other implications. 
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1 Introduction 

Gamma-ray burst (GRB) is the most explosive event 
in the universe. The standard pict ure for GRB is 



the relativistic fireball shock model (IPaczvnski 



1986 



Goodmanlll986ilShemi fc Piranlll990HRees fc Meszaros 



19921 iMeszaros fe Reeslll993t iRees fe Meszaros Ill994l) . 
Within such a picture, an initially hot fireball composed 
of photons, electron-positron pairs, and a small amount 
of baryons expands outward because of the large opti- 
cal depth, converts most of its thermal energy into the 
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bulk kinetic energy of the baryons to form a relativistic 
cold shell; the expanding shells interact with each other 
and with the surrounding medium, causing their ki- 
netic energy to be radiated in shock waves and produc- 
ing the observed GRB prompt and afterglow emissions. 
Within such a scenario, the relativistic shock gener- 
ates t he magnetic fiel d via Weibel instability (jWeible 
19591: lMedvedevlll999h and the energetic electrons via 
first order Fermi acceleration which cool down, most 



1994; 



likely via syn chrotron emission ([Meszaros et al 
Tavanilll996h . Although the standard fireball model 



can explain the general features of GRB: the early- 
time rapid temporal variability and late-time smooth 
afterglow, there are some observational features be- 
yond the expectations of this model such as the early 
X-ray p lateaus, vari ous rebrightcnings and chromatic 
breaks (jZhand 120111 ). Different extensions of the ba- 
sic model are invoked to explain the observed deviation 
from the model. The extensions include the modifica- 
tion of the total energy of the ejecta, the environment , 
the m icrophysics p arameters and the radiative mecha- 



nism (Zhang 2011). However, often these are tailored 



on a burst by burst basis. 

Most astrophysical system, e.g., accretion disks, so- 
lar/stellar winds, and the interstellar medium (ISM) are 
in turbulent states with embedded mag netic fields that 
influence almost all of their properties (Biskamp 2005 



Frischlll995HGoldstein et al 
20041 IScalo fc Elmeereen I 



1995tlelemegreen fc Scalo 
2004h . Narayan & Ku- 



mar (2009) jNarav an fc Kumar I |2009j) and Lazar et 



al. (2009) (ILazar et al. 1 120091 ) have proposed a rel- 
ativistic turbulence model instead of the internal 
shock model as the production mechanism for fast 
variable GRB light curves and ap plied it to GRB 



app 

080319B (iKumar fe Naravan 1 120091) . Zhang & Yan 



(jZhang fc Yan 11201 II ) have also developed a new model 
of GRB prompt emission in the highly magnetized 
regime, namely, the Internal-Collision-induced Mag- 



2 



netic Reconnection and Turbulence model which not 
only carries the merits, but also alleviate some draw- 
backs of the internal shock model. The role that 
the magnetic fields play in acceleration, collimation 
and emission production in GRB outflows remains one 
of the central issues. As an alternative to micro- 
scopic Weibel instability, a production mechanism of 
magnetic field by macros c opic t ur bulence is proposed 



(Go odman fc MacFadven 2007]; ISironi fc Goodman 



20071 ) and verified by recent three-dimensi onal relativis- 



tic m agnctohydrodynamics simulation (IZhang et al. 
20091) . This simulation also show that the macro- 



scopic turbulence produces small-scale moving mag- 
netic clouds which is likely a site for Fermi acceleration 
of charged particles. It is possible that the turbulence 
may give a self-consistent picture of GRB emission, 
which includes the fast variability, the acceleration of 
nonthermal electrons and the amplification of magnetic 
fields. 

The pre-shock turbulence can be triggered by the 
interaction between GRB precursors and the interstel- 
lar medium through fluid instabilities. Cosmic rays 
acceleration in relativistic collisionlcss shock can also 
excite large-scale turbulence in the shock upstream 
( Milosavlievic fc Nakar I [20061 ). In this paper, we fo- 
cus on how the turbulence/shock wave interaction mod- 
ify the propagation of the relativistic shock. Non- 
rclativistic shock wave propagation in turbulent in- 
terplanetary plasma had been studied, the result of 
which show that the solar wind plasma turbulence may 
put consid erable contribution in the shock wave de- 
celeration ( Chashei fc Shishov 1 119961 ) . In the fireball 



shock model, the interaction of ultra-relativistic ejecta 
with the surrounding medium drives a relativistic shock 
which satisfies the Blandford-McKee (BM) self-similar 
solution and totally de termines the time evolutio n of 
the afterglow emission ( Blandford fc McKee Ill976 ). If 
the turbulence exists in the upstream medium, we show 
that the relativistic shock will also amplify the turbu- 
lence behind the moving shock font, and thus deceler- 
ates more quickly than that BM solution predicted. We 
derive below the evolution solution of the relativistic 
shock propagating in a turbulent medium and discuss 
the implications for GRBs. 



2 Turbulence induced additional deceleration 
of relativistic shock 

We consider a relativistic GRB jet sweeping up the sur- 
rounding medium produced by the compact GRB pro- 
genitor with density p = Ar~ s , where A=nm p for inter- 
stellar medium (ISM) environment while A = 5 x 10 11 A* 



for wind environment ( Chevalier fc Li 2000[) . The hy- 
drodynamics involves a relativistic blast wave expand- 
ing into the wind. For an ultrarelativistic, adiabatic 
blast wave, Blandford & McKcc find that 



lirAT 2 R 3 - s c 2 
17 -4s : 



(1) 



where E is the isotropic energy of the blast wave 
(mostly denoted by E[ so ) and R is the shock wave ra- 
dius. In the observer's frame, there is t = R/(2T 2 c) 
and thus r oc i(s-3)/(8-2s)^ 

The turbulence in the wind will make T decay more 
faster. We first review the basic idea of energy transfor- 
mation process from non-relativistic shock wave to the 
turbulence, then generalize it to the relativistic shock. 
If the turbulence in the regions before and behind the 
shock surface is of acoustic type and it induced sound 
wave is incident normally from the front on a non- 
relativistic shock wave, the turbulence energy transfor- 
mation coefficient is 



X 



Wd 

w u 



2~ 



7 + 1 



Mi 



(2) 



where Wd and W u are the energy density of turbulence 
in the downstream and upstream of the shock, M u is the 
upstream Mach number, 7 is the ratio of specific heats, 
and 77 ~ 1/(7 + ^27(7 — I)). For 7 = 5/3, we have 
X — O.IM^ which means the considerable amplification 
of turbulence for strong shock waves M u ^> 1. The 
source of the turbulence amplification is the shock wave 
energy Chashei & Shishov showed that the relative 
level of turbulence Sd = Wd/pd = V 2 ^u is small, which 
means that the turbulence behi nd the shock is always 
weak dChashei fc Shishov Ill996h . 

For a relativistic shock, the Rankinc-Hugoniot rela- 
tions with the turbulence taking into account have the 
following form 



r u n u /3„ = T d n d ft d 



WuKPu + Pu + TTt, 



W d T lPd +Pd + K d 



w u TlP u + H tu = w d T 2 d /3 d + H t<i , 



(3) 



where tt u ^ are the turbulence momentum fluxes and 
-Htu.td are the turbulence energy fluxes. The rela- 
tion between itd{Htd) and n u (Ht u ) for relativistic shock 
is unclear. We just know that the energy flux can 
be approximated as H ~ Y 2 ec, where e is the en- 
ergy density of the turbulence in the comoving frame 
for relativistic flow. It is convenient to define a rela- 
tive turbulence energy transformation coefficient Sdu = 
(f/td — Ht u )/e' for relativistic shock, where e' is the 

energjj density of the post-shock medium in the comov- 
cm B 

mg frame. Due to the energy transformation from the 



3 



shock to the turbulence, the shock flow total energy 
decre ase. We should adopt an energy equation in the 
form (jChashei k Shishov I fl996h 



(4) 



where e is the energy density of the fireball (without 
the turbulence energy) in the central engine frame, H 
is the turbulence energy flux. If the turbulent terms 
are neglected in (3), the BM solution is achieved. If the 
turbulence is included, the BM solution is not valid any 
more. 

We shall use the perturbation method to solve (4) 
under the condition S d u € 1. In the central engine 
frame, we have dR = v s dt, where v s is the velocity of 
the shock wave. Carrying out the volume integration 
of the energy equation with including the energy losses 
(caused by turbulence) on the shock surface, we get 



v,— —( I 4irr 2 edr 
dR\ 



-(H td -H tu )4irR 2 



(5) 



Since J Amr 2 edr = E for zero approximation, and the 
shock is rclativistic v s ~ c, the above solution can be 
written in the form 



^ « -4n6 du e'R 2 . 
dR 



(6) 



In principle, i5d u depends on the detailed process of the 
turbulence/shock interaction. To obtain an analytical 
solution, we assume 8^ = const. 

ISM environment (s = 0). The post-shock energy 
density is e' ~ AT 2 nm p c 2 and the total energy of the 
fireball is E = 871 nm p T 2 R 3 c 2 / '17 '. Then the solution of 
(6) is 



r oc i?-3 ex p (-17 



8du 



dr 



(7) 



where i?o is the initial distance of the shock occurs. 
For the case <5d u = const, we have T oc R-z/'z-ns^ . 
Using the relation R = 2T 2 ct, the evolution solution 
of relativistic shock is T cx t-3/8-W(i6/i7+Ws) = 

£—3/8— /dec 

WIND environment (s = 2). The post-shock en- 
ergy density is e! = 4T 2 AR~ 2 c 2 and the total energy is 
E = 8ttAT 2 Rc 2 /9. We get V oc R-l/*-™* and thus 
the evolution solution is V oc t- 1 /4-5 du /(4/9+45 du ) = 

t-l/4-Atae. 

Additional deceleration factors are involved as com- 
pared with the BM solutions for ISM and wind envi- 
ronment respectively. We note that the coefficient (W 
in fact includes the energy transformation process be- 
tween the shock wave and the turbulence, which may 
be in principle defined more precisely. By defining it, 
our results can also be applicable for MHD turbulence. 



3 Implications and discussions 

In this work, we have demonstrated that the relativis- 
tic blast wave involving the turbulence in the upstream 
decelerates more quickly than what BM solution pre- 
dicted based on the condition that the relative turbu- 
lence energy transformation is small. This dynamical 
effect changes temporal behavior of the peak frequen- 
cies and thus the slope of the light curves. We re- 
strain our discussions for ISM case. The position of 
the peak of the spectrum F v varies as v c oc t~ 1 ' 2+A ^ dea 
in fast cooling while v m oc r 3 / 2_4 ^ doc in slow cooling. 
The main results are summarized in Tab. 1. Emis- 
sion for a forward shock predicts definite relation be- 
tween the spectral and t emporal indexes, F v oc t~ a v~P 



(jRees k Meszaros 1119981 ). Observed af terglow typicall y 



do not comply with this prediction ( Racusin 20091) 



The deviation of the observed closure relation from the 
external shock afterglow model in those bursts such 
as GRB050413, GRB0607A and GRB061202 may b e 
caused by the upstream turbulence (jLiang et al. Il2007h . 

The outflows of GRB are believed to be ultra- 
relativistic jets. Since the outflow decelerates more 
quickly, a jet-break in afterglow will appear earlier when 
T ~ where 9j is the jet opening angle. Using 

T oc R~3/2-ns du and T 2 R 3 = i7E iso /8Trnm p c 2 (T is 
the initial Lorentz factor of the jet), we get the jet-break 
time 



U 



1.4 



E, 



iso,53 

n 



0, ? _i(ro0j) 9+102 ^ days, 



(8) 



where 6j = O.10j ; _i, E- lso = 10 5J -Ei SO]5 3 erg and n = 
lno cm -3 . (For a GRB located at rcdshift z, the ob- 
served time should be increased by a factor 1 + z.) The 
jet-break time is earlier by a factor of the last term 
in (8) than that of the adiabatic blast wave predicted. 
For T = 1000, 6j = 0.1 and <5 du = 0.1, the value of 
the factor is about 0.2 which means that the previous 
jet opening angle is notably overestimated. The jet- 
corrected energy Ej ~ ^9 2 E iso should be potentially 
re-evaluated by (8). If the turbulence/shock interaction 
is universal, i.e., the same value of Sdu for all GRBs, 
the three parameters (tj, E iso , T ) should be statisti- 
cally correlated. Once obtain the correlation, we can 
evaluate and E^. While it ' s not easy to id entify the 

l2009h . 



jet-break (jLiang et al. I 2008; R acusin 



Because 

in most GRBs, the optical and X-ray breaks arc chro- 
matic which is contradict with the forward shock model. 
A natural but not affirmative way to remove the con- 
tradiction is to invok e a different origin of X-ray emis- 
sion: dust scattering (jShao et ah 1l2007h or photosp here 
emission of the engine activity ( Wu k Zhan j 2011 ) . So 
the late optical observation is crucial to identify the 



•1 



jet-break. Secondly, the detection of very early opti- 
cal afterglow peak can pro vide a direct m easurement 
of the initial Lorentz factor (jVergani 1120071 ). There are 
many optical afterglows without detailed X-ray cover- 
age before Swift era. Now, that Swift-XRT provides 
impressive X-ray light curves, there are too few opti- 
cal light curves. The Chinese-French mission SVOM 
is a multi-wavelength GRB observato ry scheduled to 
launch in 2014-2015 rtPaul et alj|201lh . Its operation 
window overlapping with that of Fermi and Swift will 
shade light on this problem. 

The turbulence/shock interaction make some shock 
energy stored in turbulent state, which lead to some 
uncertainties and differences between the observational 
and real values of tj, E- mo , Lj so (isotropic luminosity) 
and Speak (peak energy of vF v spectrum). A num- 
ber of relations i nvolving these parameters were pro- 
posed (see, e.g., (IFrail et al. Il200ll; lAmati et al.ll2002t 
Ghirlanda et al.l2004HYonetoku et alj2004lLiang & Zhang | 



tional Natural Science Foundation of China (grant 
11003014/A0303). 



20051) b which lead to identification of a "GRB standard 



candel" for cosmology information complementing that 
derived from SNe. This requires that the energy and 
the luminosity are precisely estimated from observa- 
tional quantities. But the significant dispersions of 
these correlations prevent GRB as a good standard 
candle, which may be intrinsic due to the lack of the 
full knowledge of the turbulence/shock interaction. 

Those quantities related to the shock will be effected 
by the turbulence, so do the correlations between these 
quantities. Although current understanding of the tur- 
bulence/shock interaction is limited, posing a challenge 
to accurate prediction of this highly non-linear phe- 
nomenon, the results show the potential ability of the 
turbulence to solve some problems in GRBs. The re- 
sults of this letter are also applicable for SNe, AGN and 
microquasar. 

Xue-Wen Liu thank Jin-song Zhao for valuable 
discussions. This work is supported by the Na- 



Table 1 Temporal index a and spectral index f3 in after- 
glow for ISM case, the convention F v oc t~ ~ a v-P is adopted. 
The power law index of the electron distribution is p > 2. 
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